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Abstract 

An approximate many-body theory incorporating two-body correlations has 
been employed to calculate low-lying collective multipole frequencies in a Bose- 
Einstein condensate containing A bosons, for different values of the interaction pa- 
rameter A = 4^ • Significant difference from the variational estimate of the Gross- 
Pitaevskii equation has been found near the collapse region. This is attributed to 
two-body correlations and finite range attraction of the realistic interatomic inter- 
action. A large deviation from the hydrodynamic model is also seen for the second 
monopole breathing mode and the quadrupole mode for large positive A. 



Low-lying collective excitations of a Bose-Einstein condensate can provide valu- 
able information about the interactions and stability of the condensate, as the dimension- 
less interaction parameter A = decreases from a large positive (repulsive interaction) 
value to negative (attractive interaction) values. Starting from the Gross-Pitaevskii (GP) 
equation, and using a sum rule approach, Stringari p] obtained analytic expressions in 
terms of expectation values of kinetic energy (Eki n ) and harmonic confining potentials 
(Eho), for the monopole (ujm) and quadrupole (u>q) frequencies. In Fig.l of Ref. [1], vari- 
ational estimates J2j [3] were plotted against A, which shows a sharp change for both ujm 
and ojq as A approaches the variational critical value \f r p ' var = —0.671 [3]. Note that 
this provides an upper bound. Furthermore, numerical calculation of Ekin by the GP 
equation close to the collapse region involves large errors. Hence, it is desirable to calcu- 
late the low-lying excitation frequencies directly using a many-body approach, which can 
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handle a finite range interaction instead of the contact interaction in the GP equation. 

In the present work, we calculate low-lying collective excitation frequencies of 
a dilute Bose-Einstein condensate (BEC), using the many-body potential harmonics 
expansion (PHE) approach [I]. Only two-body correlations are relevant in the dilute 
many-body system. Hence the (ij)-Faddeev component of the many-body wave function, 
is a function of the relative separation (fij) and a global length called hyperradius (r) 
only. In the PHE method, one expands the Faddeev component (ipij), corresponding to 
the (zj)-interacting pair of the condensate containing A atoms, in the corresponding po- 
tential harmonics (PH) basis [5] (which is the subset of the full hyperspherical harmonics 
(HH) basis [5], needed for the expansion of the two-body potential V(fij)) 

<Mw) = r^£^ + ,(%)<4(r). (1) 

Here, I and m are the orbital angular momentum of the system and its projection and Qij 
represents the full set of hyperangles for the (zj)-partition. A closed analytic expression 
can be obtained for the PH, ^ 2 K+«(^j) [5]. The expansion (1) is, in general, very 
slow due to the fact that the lowest order PH is a constant and does not represent the 
strong short-range correlation of the interacting pair, arising from the very short range 
repulsion of the interatomic interaction. To enhance the expansion we introduce a short- 
range correlation function r](rij), in analogy with atomic systems [3 [8]. This is obtained 
as the zero energy solution of the two-body Schrodinger equation with the chosen two- 
body potential [9], corresponding to the appropriate s-wave scattering length (a s ). Thus 
we replace expansion (1) by [TO] 

-(3A-4) , . 

Mr ij ,r) = r ^ SjfF^^jjXlr))]^). (2) 

Substitution of this expansion in the many-body Schrodinger equation and projection 
on a particular PH gives rise to a system of coupled differential equation (CDE) in r 
[UdO], which is solved numerically using hyperspherical adiabatic approximation (HAA) 
Details of the procedure can be found in ref. [TO] . 
In the HAA, the coupling potential matrix together with the diagonal hyper- 
centrifugal repulsion is diagonalised to get the effective potential, u; D (r), as the lowest 
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eigenvalue of the matrix for a particular value of r. Collective motion of the condensate 
in the hyperradial space takes place in the effective potential u (r). Ground state in 
this well gives the ground state energy (Eqo) of the condensate corresponding to n = 0, 
I = 0. Here E n \ is the energy in oscillator units (o.u.) of the n th radial excitation of 
the I th surface mode. Hyperradial excitations corresponding to the breathing mode for 
I = 0, give the monopole frequencies as UMn = (E n o — E 00 ). For I ^ 0, we get the surface 
modes. Numerical calculation of the off-diagonal potential matrix elements for / ^ is 
fraught with large inaccuracies and its numerical computation is very slow. On the other 
hand the diagonal hypercentrifugal term is very large for large A and contributes most to 
the potential matrix. Hence we disregard I > contributions to the off-diagonal matrix 
elements. For non-zero orbital angular momentum, we thus get the effective potential 
u>i(r) in the hyperradial space. The ground state in this potential is the ground state of 
the I th surface mode Eqi and n th radial excitations provide E n \. The lowest monopole and 
quadrupole frequencies are given by ujm = E w — Eo Q and ujq = E 02 — E 00 . The s-wave 
scattering length (a s ) has been chosen as 2.09783 x 10~ 4 o.u. for repulsive interaction and 
—1.39217 x 10 -4 o.u. for attractive interaction. The interatomic interaction is chosen as 
the van der Waals potential with a hard core of radius r c , viz., Vfaj) = oo for < r c 
and = — et for > r c . The value of C§ is chosen to be appropriate for rubidium atoms. 
The value of r c is adjusted to get the desired value of a s [12]. All quantities are expressed 
in oscillator units (o.u.) appropriate for the JILA experiment with 85 Rb [13J. The PHE 
equation is then solved numerically for different A to calculate E n i for different values of 
A. 

In Fig.l we plot calculated ujm and uoq as a function of the dimensionless inter- 
action parameter A. For comparison we include the variational estimates from the GP 
equation [lj. We have also included results of a numerical solution of the GP equation 
[T3] , However calculation of Eki n and Ei nt using the GP program have large errors for 
values of A close to collapse for attractive condensates as also for large positive A. This 
is shown by the lack of agreement with the virial identity involving E^ w E int and Eh , 
viz. 2Ekin — 2Eho + SEint = [15J. Hence in Fig.l, the curves obtained by a numerical 
solution of the GP equation is shown only in the region where it is reliable. It is seen that 
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for positive values of A, all the three curves for the monopole frequency are close to each 
other. Results of PHE and variational estimates are almost indistinguishable. However 
for negative A, uom by the PHE method drops faster as A approaches the critical value 
A cr = —0.461. Note that while GP equation predicts a critical value of —0.571 (varia- 
tional estimate is —0.671) [15], experimental critical number is — 0.459±0.012±0.054 [13j. 
The prediction by the PHE method agrees well with this value [16]. As the attractive 
condensate approaches collapse, with A approaching A cr from the right, the condensate 
becomes softer (more compressible) and the breathing mode easier. 
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Figure 1: Calculated monopole (w«) and quadrupole (uq) frequencies as functions of A. 
Superscripts GP, GP,var and PHE correspond respectively to numerical solution of the 
GP equation, its variational estimate, and the present calculation. 

On the other hand the quadrupole frequency {ujq) shows a greater deviation from 
the variational estimate - both for positive and negative values of A. Since A cr = —0.461, 
it is natural that both monopole and quadrupole frequencies calculated by PHE method 
have rapid changes as A approaches —0.461. The difference between the collective excita- 
tion frequencies calculated by PHE and GP equation can be attributed to the two-body 
correlations and finite range interactions. As the condensate approaches collapse, its 

4 



size decreases rapidly |16j . so that the finite range interatomic attractions of the realistic 
atom-atom potential dominate and correlations play a more important role. Although for 
small values of |A| use of an effective interaction in terms of a s is adequate, the strong cor- 
relations arising from the attractive part of the realistic atom-atom interaction produces 
discernible effects [TTJ HJ [191 [20] f° r ^ approaching A cr . One can also qualitatively inter- 
pret the dependence of ujq on A as follows. The two-body correlation disallows proximity 
of two interacting atoms in the repulsive case. Again in the I = 2 state the atoms are 
further apart compared to the / = state. It is therefore easier {i.e. excitation energy is 
lower) for the condensate to jump from / = to / = 2 state, when two-body correlations 
are included than when they are excluded (as in the GP equation). The inverse is true 
for the attractive case. The results obtained in the present work are in conformity with 
these arguments. Further discussions on the behaviour of ujq for large positive A will be 
taken up together with those of higher excitation frequencies, after we present them. 
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Figure 2: Higher order multipole frequencies calculated by PHE method, as functions of 
A. 

In Fig. 2 we present the multipole frequencies which are degenerate with a 
value of 4 in the non-interacting case, viz., u> M2 = E 2i o — E , u>q 2 = Ei )2 — E Q and 
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ujh = E 0} 4 —Eo,q- As expected from the softness of the condensate , the second breathing 
mode frequency ujm2 is larger than 4 for positive A and less than 4 for negative A, de- 
creasing sharply as A — > X cr . This once again shows that the condensate becomes softer 
as A — > A cr , and harder and more rigid as A increases to large positive values. On the 
ther hand, both ujQ2 and ujh are less than 4 for A > and greater than 4 for A < 0. 
Both increase sharply as A — > X cr , ujh increasing much faster than ujq2- As positive 
A increases, the condensate becomes less compressible (more rigid) costing less energy 
for rotational excitations. Hence the spacing between the rotational levels decrease. For 
negative A close to A cr , uh is much larger than uq2- This shows that for such soft conden- 
sates, the first breathing mode of I = 2 surface mode is easier to excite than the ground 
state of the I = 4 surface mode, although they are degenerate in the non-interacting 
limit. This is consistent with the general picture and intuition. One can also note that 
uq2 ~ = Ei t 2 — E 0) 2 decreases rapidly from 2 as A decreases from zero to A cr . This 
again shows the softness of the / = 2 surface mode of the attractive condensate. 

We next come back to a discussion of our results as compared with those by the 
GP equation and the hydrodynamic model, as also with available experimental results. 
As mentioned earlier, results using the numerical solution [T4] of the GP equation are 
not reliable for almost all values of A except small |A|. Moreover ujm and ujq given in 
terms of E\~ in and Eho are rigorous upper bounds for the frequency of the lowest states ex- 
cited by the appropriate multipole operator pQ. The upper bounds get further upwardly 
shifted for the variational estimate. The PHE result for ujq is lower than the variational 
estimate; hence it is consistent with this requirement. Using the hydrodynamic model 
Stringari also derived a dispersion law for the frequencies of the normal modes p], viz., 
u(n, I) = co> (2n 2 + 2nl + 3n + 1) 1 / 2 , where Uq is the frequency of the spherically symmetric 
trap. This is valid in the limit A >> 1, where kinetic energy contribution can be disre- 
garded. While the PHE result for ojm is very close to this estimate, that for ujq is close 
to but less than this limit. Available experimental results [21| [22| [23j [23] are for axially 
symmetric traps only - measured uj for m = ±2 and m = agree well with the GP 
equation and the hydrodynamic model estimates for the deformed trap. Since accuracy 
of the hydrodynamic model estimate is expected to decrease as n and / increase, due 

6 



to the neglect of the kinetic energy contribution pQ and ujm calculated by PHE agrees 
remarkably well with that from the GP equation, it is important to have accurate mea- 
surement of ujq for the spherical trap. The experimental procedure of ref. [21] can be 
utilized for this purpose, since it is suitable for any geometry of the trap. This technique 
uses the optical dipole force of a rapidly scanning laser beam to excite surface modes by 
inducing deformations of the trap potential. 

More remarkable is the disagreement of the second breathing mode (umz)- The 
PHE predicts a value > 4 for positive A, while the A >> 1 limit by the hydrodynamic 
model gives y/Ti. One expects that the system becomes less and less compressible as 
positive A increases. Consequently one would expect um2 > 4 for positive A, which dis- 
agrees with the hydrodynamic model prediction. This is an indication that the kinetic 
energy contribution can no longer be disregarded for such excitation energies. To resolve 
the issue, one needs a precise measurement of a>M2- The higher mode can be excited by 
an appropriate frequency of the laser beam. 

In summary, we have investigated monopole, quadrupole and hexadecapole fre- 
quencies (corresponding to I = 0, 2 and 4 respectively) of the collective motion of 
Bose-Einstein condensates by a many-body approach based on hyperspherical harmonics 
method. We include the most important two-body correlations, disregarding higher-body 
correlations. We notice that two-body correlations and the finite range attraction of a 
realistic two-body interaction make the attractive condensate softer compared with the 
use of a contact interaction. This is consistent with intuitive expectations. Consequently 
our results differ appreciably from those by the GP equation in the collapse region for 
an attractive condensate. The effect of two-body correlations become increasingly more 
relevant as the condensate approaches collapse. For repulsive condensates the monopole 
frequency agrees very well with the the GP estimate, while the quadrupole frequency 
differs by about 10% of the asymptotic value. For the second breathing mode the PHE 
result is greater than the non-interacting value, as expected intuitively, for A > 0. This 
is in sharp contrast with the asymptotic (A >> 1) GP estimate of V44, indicating that 
the disregard of kinetic energy terms in the hydrodynamical model is not valid for such 
excitation levels. This calls for precise measurements of these multipole frequencies. 
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